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Abstract 

By describing a large class of deep inelastic processes with standard parameterization for the dif- 
ferent parton species, we check the characteristic relationship dictated by Pauli principle: broader 
shapes for higher first moments. Indeed, the ratios between the second and the first moment and 
the one between the third and the second moment for the valence partons is an increasing function 
of the first moment and agrees quantitatively with the values found with Fermi-Dirac distributions. 



Four experimental facts conspire to indicate that the Pauli exclusion principle plays a role in the 
quark-parton distributions in the nucleons. They are the defect |l| in the Gottfried sum rule ||], the 
asymmetry in the Drell-Yan processes on proton and deuteron targets |J, and the high x behaviour 
of the ratios F^ix^/F^ix) [f| and cfi(x)/Ff(x) ||. The first two facts imply the inequality d > u 
in the sea of the proton, which has been advocated long time ago [||. From the other two one can 
deduce the dominance of the u* parton in the high x region; this is an indication of the relationship 
between the shape of the distribution as a function of x and the first moment of a parton required by 
Pauli principle: broader shapes for higher first moments. 

More recently, by assuming the approximate relationship u^(x) = —d(x), it has been possible to 
successfully relate the structure functions F^{x) — F^{x) and xg\(x) for x > 0.2. 

Pauli principle suggests to assume [|| Fermi-Dirac functions for the quark distributions, 

exp{ 1 '-} + 1 

and Bose-Einstein functions for the gluons, 

°exp{ ^ '-} - 1 

where x plays the role of the temperature, the x, are the thermodynamical potentials of each parton 
species and f(x) is a weight function with the usual form A x a (1 — x) 13 . 

To account for the increase of F%{x) at small x |J we have to add to Eq. (jl|) an additional 
contribution Ql(x) = Ql(x), which should be unpolarized and isospin- invariant to get finite quark- 
parton model sum rules (QPMSR). This leads to a satisfactory description of a large class of deep 
inelastic data in Ref. [jl0| and its updated version [O], including the E154 data by SLAC [12]. 



Here we want to test a specific property of the quantum statistical distributions given by Eq. (|T]), 
namely the shape-first moment relationship previously mentioned for fermionic partons. The QPMSR 
imply 

u T » d l > u l >(i T . (3) 
The defect in the Gottfried sum rule yields 

d > u . (4) 

We can argue the same pattern than Eq. (^) and ([|) for the corresponding broadnesses. To be more 
precise we expect the ratio between the second and the first moment, that is the mean value of x, as 
well as the one between the third and the second moment of each quark parton to be an increasing 
function of the first moment. For the gluons the deviations from the Boltzmann distribution (implying 
the shape to be independent from the first moment) are in the opposite direction; thus, one expects 
AG(x) softer than G(x), and G(x) softer than q(x), if the same f(x) is taken for quarks and gluons. 
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To establish whether experiments confirm this shape-first moment relationship, we try to describe 
the same data considered in , by taking, for each light quark-parton species of a given flavour and 
spin, at Qq = 3GeV 2 , the form (p = u* , u^, d* , d^, v) = u*, d) = d^)\ 

p(x, Ql) = A p x a (l - xp + A L x a t (1 - . (5) 

Going on with the thermodynamical language, we name gas and liquid the two terms in Eq. (||). 

We are not very sensitive to the strange quark distribution. Indeed, the only observables depending 
on it that we are considering are FQix^/F® (x) and F^ix). For the first one isospin invariance implies 
the same contribution for both members of the ratio, the latter being near to one in the small x region, 
where strange partons are concentrated. In that region the second observable, F^ix), receives a large 
contribution from the second term of Eq. (||), responsible for its increase at small x. Therefore, in 
order to avoid the introduction of new parameters, we take, according to Ref. [13], 



u(x) + d(x) , . 

S ( x ) = s(x) = V } A 2 K 1 . (6) 
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As in Ref. |l| we consider F^(x), measured at Q 2 = 3GeV 2 Q, F%(x) and F 2 n (x) at Q 
AGeV 2 , |l|, ^, [0| and the polarized structure functions measured at SLAC {g\{x) and <7i(x) at 
Q 2 = 3GeV 2 |R|, gf(x) at Q 2 = 2GeV 2 and Q 2 = 5GeV 2 @). To report all the data at the 
same Qq = 3GeV 2 we assume the Q 2 dependence found by NMC, while the good agreement, in the 
common x range, between the data regarding g™(x) at Q 2 = 2GeV 2 and at Q 2 = 5GeV 2 makes us 
confident that they both give a good approximation of that quantity at Q 2 = 3GeV 2 . 

The unpolarized structure functions are given by: 

f 1 - 5 

F% ,n (x) = x < ± - [u(x) + u(x) — d(x) — d(x)] + — [u(x) + u(x) + d(x) + d(x)] 
t 6 18 

+ I [ 8 (x)+s{z)]\ , (7) 
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x F3(x) = x [u(x) + d{x) + s(x) — u(x) — d{x) — s{x)] . 
As long as for the polarized distribution we recall that the presence of a gluon isoscalar contribution, 



related to the anomaly, has been advocated [18] to explain the defect in the Ellis- Jaffe sum rule [19] 
for the proton found by the EMC experiment ||. This contribution, up to first order in a s , is given 
by the convolution |2(| 

In the framework of a complete NLO analysis, we should consider all possible a s contributions, but this 
is beyond the purposes of the present work. However, in a factorization scheme where chiral symmetry 



is respected (a choice closer to the picture of the nucleons in term of constituent quarks [21Q, the high 
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value of AG required to explain the experimental results suggests that the gluon polarization gives 
the dominant contribution proportional to a s . Therefore, we shall give the following expressions for 
the nucleon polarized structure functions: 

<??'» = ± 1 [Au(x) - Ad(x)] + A [Au (x) + Ad(x)] + S(x) , (9) 

with 

AG(x, Ql) = A AG x a ^ (1 - xf AG . (10) 

Eqs. (||) and ( ]l0| ) give a three parameter description of S(x) (a s A A g, aA& an d Pag), which is 
general enough for the gluon isoscalar term that one has to add to the valence term to obtain g\ (x) 
and gx(x). However, the fact that the presence of this contribution, which is important in the small 
x region, influence only a little the determination of the shapes of the valence quarks, which depend 
on the /3p's, makes the debate on this issue purely academic with respect to the main purpose of this 
paper. 

As usual, we assume 

gi(x)= 1 -(l-^ D y g p 1 (x)+g n l (x)), (11) 

where lod is the amount of D-wave in the deuteron ground state (0.058) | J22| . 

We have then four parton distributions (it, d, d and u) to describe the three unpolarized structure 
functions (F% (x), F£(x), and F%(x)), and three parton distributions (Aw, Ad and AG) for the two 
polarized ones (g^ix) and g™(x)). To overcome the ambiguity arising from this situation, we require the 
unpolarized distributions to comply with the Adler sum rule H| and with the NA51 g measurement 
of the asymmetry in the Drell-Yan production of muon pairs, 

a DY (x = .18) = -0.09 ±0.04 ±0.02. (12) 

We also require the (5q and Pag to be larger than all the P q : s and 

- 0.8 < a, a A G < . (13) 

We consider the cases with AG = and with AG free. At small x it is difficult to separate the two 
terms in Eq. (||). Since the second term is expected to give a small contribution (~ 10%) [|l(]] to 
the second moment of each parton, but an infinite one to the first moment, we expect the second 
moment of the first term to be less affected than the second one by the ambiguity in disentangling 
the two contributions. However, the second term gives an equal contribution to each quark-parton: 
therefore the existence of a correlation between the contribution to the first moment and the shape of 
the first term in Eq. (||) is not affected by this uncertainty, which has the same influence on the first 
and second moment of the gas component of each parton. Moreover, since the low x behaviour of the 
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two contributions is very different (in the case of Ref. [1C] their power behaviours differ by almost 
one unity), we expect that also for the first moment the ambiguity is not so large. 

As we have previously said, the Fermi-Dirac function for the quark distributions gives an automatic 
correlation between the first moment and the shape of each quark-parton distribution, namely broader 
shapes for higher first moments, with both quantities depending on x, once fixed f(x) and x. Instead, 
for the first term of Eq. (|5|), once imposed a common value for a, there are two parameters for each 
parton, A p and (3 p , to reproduce the first two moments of the first term in Eq. (||), which we shall 
indicate with p and p^ respectively. The second term in Eq. (||), responsible of the increase of F%{x) 
at low x and apparently power-like, gives an infinite contribution to any p (equal for each light parton, 
should we start from a finite positive value of x) and an equal contribution to p( 2 \ 

Since it is difficult to disentangle the gas and liquid component in Eq. (||), we parameterize each 

light parton distribution with the same a; in this way the shape of the corresponding distribution is 

dictated by (3 P , which is related to the high x behaviour, where the contribution of the liquid part is 

expected to become negligible. To make the comparison between Eqs. (jlj) and (|5|) more appropriate, 

we fix the liquid part to the value previously found with the parton statistical distributions, namely 

0.12 x- 1 - 19 (1 -x) 9 - 8 AG = 0, 
A L x aL (1 - xf L = (14) 
0.12 x~ L19 (1 -x) 9 - 6 AG free. 

In tables ffl, ^ and ^ the values of the parameters found in the fit and the first three moments 
of the gas component, with their ratio for each parton, with and without gluons, are reported and 



compared with the updated fit (whose results are almost identical to the one in Ref. [11]) with 
quantum statistical distributions. 

An attentive glance to the numbers in the tables brings us to the following conclusions: 

i) for the quantum statistical as well as for the present parameterization the parameters and the 
moments of the quark-partons have a very weak dependence on the gluon contribution and the 
agreement between the moments, obtained with the different parameterization, is good with the 
exception of the u parton and AG. 

ii) using parameterization (||) we get, for the ratios pC 2 ' /p and p( 3 )/p( 2 ) as a function of p for the 
quarks, a pattern very similar to the one obtained with Fermi-Dirac functions, which have the 
property 

«T(2) rfi(2) u m rfT(2) 

— » — ^ — >^T' ( 15 ) 

U T(3) d im > n J(3) d T(3) 



» "TTToT" £ TToT > 



„T(2) " rfi(2) ~ u i(2) " dT( 2 ) ' 

this shows that the ratios p^ jp and p( 3 )/p( 2 ) are increasing functions of p according to Eq. (]3|). 
Also, S 2 ^/d> vP> /u, as we expect from Eq. @. 
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The differences found for u and AG are not unexpected. As for the first one, which is the lighter 
parton with the lowest first moment, the difference is large in percentage but small in amount. The 
different shape found, which is a consequence of the high value of /3a, is explained by the small 
precision in the determination of a parameter which is related to the shape of a distribution with a 
very low first moment. Concerning the gluons, we are not surprised to find a AG about more than 
twice larger than with the Bose-Einstein form, which has a less divergent extrapolation to small x, 
and a different shape: the fact that the fit is equally good with and without gluons shows that we 
are not able to get, within our approach, relevant information on them. However, the value found is 
consistent with the NLO result of Ref. p5| , AG = 1.1 ± 0.4, where the gluons are also constrained by 
the Q 2 dependence of the data. The value found for ckag is & t the lower limit of the one in Ref. [25]. 

The fact that the ratios S^ 2 ' /S and d)^ /cft^ 2 ', for the valence quark with the smallest first 
moment, are smaller than the ones for the other valence quarks confirms the correlation among 
broader shape and higher first moment, which shows up in the dominance of v} at high x. Needless 
to say, iJ( 2 ) /v} and u)^ /v}^ come out larger than the corresponding ratio of the other valence 
partons, as we expect from the high x behaviour of F^x) / F%{x) and A±(x). Indeed, already in 



Ref. [26], by taking the single parton distributions from a set of deep inelastic data, the author has 
been able to find a similar trend with a distribution for v) and for & broader and narrower respectively 
than for the other valence partons. 

The values of the first and second moment of the gas part in Eqs. (Q) and (|5|) are in good agreement. 
(The exception of u does not contradict this conclusion, since the small value of its first moment makes 
less relevant its shape.) This agreement is, at our advice, an additional argument in favour of the 
Fermi-Dirac distributions, for which, once fixed x and f(x) for each parton, p, p^ 2 \ and p^> depend 
on one parameter, x(p); to say it in a more clear way, the fit with the present parameterization might 
give equal values of /3„'s and distributions which differ only by a constant factor, which would be, 
within our statistical language, the Boltzmann limit. The agreement found for the ratio p^ /p^ for 
the valence quarks with the two different parameterizations gives a positive test, which is very little 
affected by the difficulty of disentangling the gas from the liquid contribution, since the last one has 
little relevance in the evaluation of the second moments of the valence quarks and even a minor one 
for the third moments. 

We wish to compare our evaluation of the first two moments of the non-singlet polarized distribu- 
tions with a previous analysis of the same data [27] and with the result found in the NLO analysis of 
Ref. [|25[|. To this purpose, in table || these moments and their ratio are reported for g^ N . One can 
see a very good agreement of the value found for the second moment of g±^ NS ^ between the quantum 
statistical distributions, the present parameterization and the one in [p5f |, while the difference found 
for the first moment depends on the small x behaviour. 
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The gluons are expected to have a softer distribution than the <f s if the same f(x) is taken. The 
great uncertainty on their distribution does not allow to draw meaningful conclusions from the fact 
that table || agrees with this expectation. A typical property of the Bose-Einstein function is that 
AG(x)/G(x) is a decreasing function of x (soft polarization for the gluons), therefore AG^/AG < 
G^ /G. In order to show this property one needs a NLO analysis suitable to determine the shape of 
AG and G. We plan to perform a consistent NLO analysis to include all the contributions proportional 
to a s (Q 2 ) and to evolve the distributions found in this way to the higher Q 2 (~ 10 GeV 2 ) of the CERN 
experiments on polarized proton and deuteron targets. We do not expect a too challenging test for 
the resulting predictions, since the precision of the measurements, especially in the low x region not 
reached at SLAC, is worse for CERN data. 

In fact, we tried to extend our analysis to the CERN data at <Q 2 >= 10 GeV 2 for g\(x) and 
gf(x) p8|| , together with the unpolarized distributions evaluated at Q 2 = 10 GeV 2 . Unfortunately 
the absence of data on g™ (g^ 3 ) and the minor precision of the data for gf, together with the fact 
that the unpolarized structure functions at Q 2 = 10 GeV 2 are known with a smaller precision, do not 
allow to extract from the data the trend of p^ jp and p^ /p^ , which come out different from equally 
good fits. 

To get an idea of quality and differences of the two fits we compare in Fig. |l| and |2| x F^{x) 
measured in Ref. [|14|] and x gi(x) measured in Rcf. [12] with the quantum statistical and present 
parameterization in the case AG free. As we see, the data are very well reproduced and there is very 
small difference between the two fits. 

In conclusion, as a consequence of our analysis, we can add a new fact to the ones mentioned at 
the beginning of this paper, supporting a role of the Pauli principle: S , the valence quark with the 
lowest first moment, is also the one with the narrower shape, and the dependence of the shape on the 
first moment for each parton is in good agreement with the one found using Fermi-Dirac functions. 

Our preliminary work on NLO analysis of the SLAC data confirms the trend with a slightly 
narrower distribution for d) and, of course, a largely broader one for u* with respect to and 

The study of the Q 2 evolution is in our plans, especially since the large Q 2 and x excess found at 
Hera p9[ ] might be explained by a modification of the evolution, which is resonable if Pauli principle 
plays a role in parton distributions, consisting in a slower narrowing of the distribution of a fermionic 
parton with its levels almost completely occupied (in the case of the proton the v}), as advocated in 
Ref. [001. 
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TABLE 





PP 




ST 




AG = 


AG free 




AG = 


AG free 


a 


-0.227 ±0.012 


-0.235 ±0.012 


A 

a 

P 


o ak + n f)7 

i-j . l 'J _i_ yj.KJi 

-0.252 ±0.009 
2.19 ±0.04 


o so + n 07 

£j.%J\J _1_ u.U i 

-0.251 ±0.005 
2.20 ±0.04 


PuT 


2.35 ±0.04 


2.35 ±0.04 




1 


1 


Pdl 


4.21 ±0.10 


4.12 ±0.11 


x(d l ) 


0.206 ±0.006 


0.182 ±0.011 


Pui 


4.35 ±0.16 


4.26 ±0.17 


x(u^) 


0.134 ±0.009 


0.112 ±0.013 


P d r 


4.81 ±0.24 


4.89 ±0.25 


x(d^) 


-0.072 ±0.012 


-0.073 ±0.013 




4.81 ±0.10 


4.90 ±0.10 


x(d) 


-0.351 ±0.016 


-0.388 ±0.024 


Pu 




i6-s±?:I 


x(u) 


-0.647 ±0.031 


-0.715 ±0.046 


CtAG 




-0.80 ±0.61 


x(G l ) 


-0.129 ±0.009 
-0.129 ±0.009 


-0.119 ±0.010 
-0.191 ±0.029 


Pag 




20 ±11 


X 


0.230 ± 0.006 


0.243 ±0.006 


X 2 


1.89 


1.90 


2 

X 


2.42 


2.37 


Bj 


0.163 ±0.003 


0.162 ±0.003 


Bj 


0.165 ±0.010 


0.164 ±0.009 


Gott 


0.223 ±0.006 


0.223 ±0.006 


Gott 


0.209 ±0.012 


0.209 ±0.015 



Table 1: The values of the parameters found with the present parameterization (PP) are compared 
with the ones obtained using quantum statistical distributions (ST). The liquid part is given by 
Eq. 0. In the last two rows we report the results for the rhs of the Bjorken |24j and Gottfried 
sum rules. 
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TABLE | 



AG = 


PP 




p 


p (2) 


p(2) jp 


p (3) 


p(3) /p(2) 




1.246 ±0.010 


0.234 ± 0.004 


0.187 ±0.003 


0.081 ± 0.002 


0.346 ±0.010 


d l 


0.664 ±0.010 


0.086 ± 0.002 


0.129 ±0.004 


0.0218 ±0.0009 


0.254 ±0.012 




0.593 ± 0.009 


0.075 ± 0.002 


0.126 ±0.005 


0.0186 ±0.0010 


0.25 ±0.02 


s 


0.342 ± 0.009 


0.040 ± 0.002 


0.117 ±0.006 


0.0094 ± 0.0007 


0.23 ±0.02 


J/2 


0.137 ±0.004 


0.0160 ± 0.0006 


0.117 ±0.005 


0.0037 ± 0.0002 


0.234 ±0.013 


u/2 


0.054 ± 0.005 


0.0023 ± 0.0003 


0.042 ± 0.007 


0.00021 ± 0.00005 


0.09 ±0.02 


ST 




P 


p (2) 


p(2) jp 


p(3) 


p(3) /p(2) 




1.25 ±0.05 


0.228 ± 0.008 


0.183 ±0.009 


0.078 ± 0.003 


0.34 ±0.02 




0.68 ±0.03 


0.091 ± 0.003 


0.134 ±0.007 


0.0240 ±0.0010 


0.26 ±0.02 




0.59 ±0.03 


0.076 ± 0.003 


0.128 ±0.008 


0.0194 ±0.0010 


0.26 ±0.02 




0.35 ±0.02 


0.041 ± 0.002 


0.116 ±0.010 


0.0097 ± 0.0007 


0.24 ±0.02 


J/2 


0.134 ±0.011 


0.0144 ± 0.0014 


0.107 ±0.014 


0.0033 ± 0.0003 


0.23 ±0.03 


u/2 


0.041 ± 0.006 


0.0043 ± 0.0007 


0.10 ±0.02 


0.0010 ± 0.0002 


0.22 ±0.05 



Table 2: The values of the first, second and third moment and of the ratios /p and p^ /p^ for 
the gas component in Eq. (||) are reported in the case AG = 0. 
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TABLE | 



AG free 


PP 




P 


p (2) 


p( 2 ) lp 


p (3) 


p(3) / p (2) 


«t 


1.254 ±0.010 


0.233 ± 0.004 


0.186 ±0.004 


0.080 ± 0.002 


0.345 ± 0.010 


d l 


0.661 ±0.011 


0.086 ± 0.002 


0.130 ±0.004 


0.0220 ± 0.0009 


0.256 ±0.013 




0.592 ± 0.009 


0.075 ± 0.003 


0.127 ±0.005 


0.0189 ±0.0011 


0.25 ± 0.02 


s 


0.350 ±0.011 


0.040 ± 0.002 


0.115 ±0.007 


0.0093 ± 0.0007 


0.23 ±0.02 


J/2 


0.138 ±0.004 


0.0158 ±0.0006 


0.115 ±0.005 


0.0036 ± 0.0002 


0.230 ± 0.013 


u/2 


0.055 ± 0.005 


0.0023 ± 0.0003 


0.041 ± 0.007 


0.00021 ± 0.00005 


0.09 ±0.02 


AG 


1.9 ± 1.3 


0.018 ±0.056 


0.009 ± 0.030 


0.0010 ± 0.0036 


0.05 ±0.26 


ST 






P 


p (2) 


p( 2 )/p 


p (3) 


p (3) / p (2) 




1.26 ±0.04 


0.229 ±0.008 


0.182 ±0.009 


0.078 ± 0.003 


0.34 ± 0.02 


d l 


0.66 ± 0.02 


0.088 ± 0.004 


0.134 ±0.008 


0.0235 ± 0.0011 


0.27 ±0.02 




0.57 ±0.02 


0.074 ± 0.004 


0.129 ±0.008 


0.0193 ±0.0011 


0.26 ± 0.02 


£ 


0.36 ± 0.02 


0.044 ± 0.003 


0.119 ±0.009 


0.0107 ±0.0007 


0.25 ± 0.02 


J/2 


0.131 ±0.013 


0.014 ±0.002 


0.11 ±0.02 


0.0034 ± 0.0004 


0.23 ±0.04 


u/2 


0.038 ± 0.008 


0.0040 ± 0.0008 


0.11 ±0.03 


0.0009 ± 0.0002 


0.23 ±0.07 


AG 


0.91 ±0.04 


0.063 ± 0.003 


0.070 ± 0.004 


0.0113 ± 0.0006 


0.179 ±0.012 



Table 3: Same as table ^ with AG free. 
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TABLE | 





PP 


ST 


Ref. H 


Ref. f§ 


p 


0.0908 


0.0927 


0.0996 


0.0989 




0.0201 


0.0197 


0.0201 


0.0226 


p(2)/ p 


0.221 


0.213 


0.202 


0.228 



Table 4: The values of the first and second moment of the non-singlet part of with AG free, ob- 
tained with the present parameterization (PP), are compared with the same quantities from quantum 
statistical distributions (ST) and from Ref. |g5[ |27fl. The values from Ref. [25] have been reported to 
Q 2 = 3 GeV 2 with the NLO evolution equation for the moments. 
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Figure 1: The prediction for xFs(x) at Q 2 = 3 GeV 2 in the case AG free is plotted and compared 
with the experimental data |14|] . The solid and dashed line corresponds to the fit with present 
parameterization and with quantum statistical distribution respectively. 
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Figure 2: The prediction for xg™(x) at Q 2 = 3 GeV 2 in the case AG free is plotted and compared 
with the experimental data 113]. The solid and dashed line corresponds to the fit with present 
parameterization and with quantum statistical distribution respectively. 
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